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1 Introduction 

Many applications use processes that are described by stochastic differential equations (SDE). 
Recently, much attention has been paid to SDEs driven by the fractional Brownian motion 
(fBm) and to the problems of statistical estimation of model parameters. Statistical aspects 
of the models including fBm were studied in many articles. Especially much attention was 
paid to the estimations of the parameter of drift. We concentrate on estimation of the Hurst 
parameter and volatility of SDE with polynomial drift 

X t =Z+ [ (aXr + bX s )ds + c f X s dBf, £>0, t € [0,T], (1) 

Jo Jo 

where 1/2 < H < 1, m £ N, m ^ 2, and B H is a fBm with Hurst index 1/2 < H < 1. Almost 
all sample paths of B H have bounded p -variation for each p > l/H on [0,T] for every T > 0 
(a more detailed definition of an fBm will be formulated in the next section). The second 
integral in 0 is the pathwise Riemann-Stieltjes integral with respect to the process having 
finite p-variation. 

It is well known that equation 0 for standard Brownian motion, i. e. for H = 1/2, is 
nonlinear reducible SDE with polynomial drift of degree m. Under the substitution h(x ) = 
x 1 ~ m equation 0 reduces to a linear SDE with multiplicative noise. Some famous equations 
can be obtained as particular cases of this equation. For example, we can obtain the Black- 
Scholes model, the Verhulst and the Landau-Ginzburg equations. 

Using an approximation approach in L 2 Dung ;7] obtained an explicit solution of the 
equation 0 for H > 1/2. His approximation is based on the fact that fBm of Liouville form 
can be uniformly approximated by semimartingales. 

Many authors (see [81, [9., 2], [5], pQ, [2], [16], [14], Q21) considered an almost sure con¬ 
vergence and an asymptotic normality of the generalized quadratic variations associated to 
the filter a (see M) of a wide class of processes with Gaussian increments. Conditions for this 
type of convergence are expressed in terms of covariances of a Gaussian process and depend 
on some parameter 7 G (0,2). These results make it possible to obtain an estimator of the 
parameter 7 and to study its asymptotic properties. If the process under consideration is the 
fBm then this parameter is the Hurst index 0 < H < 1. 

*This research was funded by a grant (No. MIP-048/2014) from the Research Council of Lithuania. 
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There are a few results concerning estimation of the Hurst index if the observed process 
is described by SDE driven by fBm. For the first time, to our knowledge, an estimator of the 
Hurst parameter of a pathwise solution of a linear SDE driven by a fBm and its asymptotical 
behaviour were considered in [4]. A more general situation was considered in D3H using a 
different approach. Estimators for the Hurst parameter were constructed making use of the 
first and second order quadratic variations of the observed values of the solution. Only the 
strong consistency of these estimators was proved. A comparison of different estimators of 
the Hurst index for linear SDE driven by fBm was carried out in El- A rate of convergence 
of these estimators with probability 1 to the true value of a parameter was obtained in m, 
when the maximum length of the interval partition tends to zero. 

In this paper we will use a different approach to solving the equation 0 than suggested 
by N. T. Dung. We apply a pathwise approach to stochastic integral equations and use a 
chain rule for the composition of a smooth function and a function of bounded p-variation 
with 1 < p < 2. This approach allows us to simplify the proof of existence and uniqueness 
of the solution of the equation (1). The goal of the paper is to construct strongly consistent 
and asymptotically normal estimators of the Hurst parameter and the volatility of SDE 0 
using discrete observations of the sample paths of the solution of the SDE. Estimators are 
constructed using quadratic variations of second order increments. Estimators’ properties are 
determined asymptotically from the properties of quadratic variations. 

The paper is organized in the following way. In Section [2] we present the main results of 
the paper. Section [3] is devoted to several known results needed for the proofs. In Section [4] 
we prove the uniqueness of SDE <[T|) and give its explicit solution. Sections 5.1 5.3 contain 
the proofs of the main results. Finally, in Section [ 6 ] some simulations are given in order to 
illustrate the obtained results. 


2 Main results 


Let X be a stochastic process satisfying 0. Denote 


where X n ^ = X(t £), = YL. For the sake of simplicity, we shall omit the index n for points 

t% and for notions X n ^ : , A^,, A„ : k, where the sample size is equal to n. In case of different 
sample sizes the indices will be retained. 

In the formulation of Theorem [2] as well as in the sequel we make use of the symbols O r ,o r . 
Here is a brief explanation. Let ( a n ) be a sequence of real numbers. Y n = O r (a n ), a n i 0, 
means that there exists a.s. finite r.v. with the property \Y n \ ^ c • a n \ Y n = o r (a n ) means 
that \Y n \ ^ c • b n with b n = o(a n ). Particulary Y n = o r ( 1) denotes the sequence (Y n ) which 
tends to 0 a.s. as n —> oo. 

Theorem 1 Assume that X is a solution of SDE 0 with a ^ 0. Let Xo > 0, c ^ 0, and 
1/2 < H < 1. Then 


2 ^ In £ ( H « 


—¥ H a.s., 
-H) Af(0,cr 2 ) 


( 2 ) 


with a known variance a 1 = cr“{H) defined in Subsection 3.2 where for n > T 


H, 




Aj 2) A- 

cXk 


ip n is the inverse function of <p n ,T, x € ( 0 , 1 ). 

Theorem |T| gives an estimator of H under the assumption that c is known. This is not always 
the case. Therefore we present another estimator which is suitable when c is unknown. 




and 
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Theorem 2 Let assumptions of Theorem^ hold and 



Then 


II 


( 2 ) 


H a.s., 


2 In 2 Vn(^n 2) — H) Af(0, a 2 ) 

with a known variance cr 2 = al(H) defined in Subsection 


( 3 ) 


3.2 


To estimate c one can use the last theorem provided below. 

Theorem 3 Assume that 1/2 < H < 1 and = H + o r (</(n)) a.s. If fi(n) = o then 

>(3)_ —1 / . fO\ \ 2 


n 2H^-l g f A ^) X 


T 2Hi?\ 4_ 2 2iri 3) ) ^ 


x k 


a.s.; 


if 4>{n) = o ^ v<g 1 lnn ) then sfhlfil - c 2 ) A/"(0, c 4 ct 2 ). 

Hence estimation of c and construction of confidence intervals is always possible. However 
confidence intervals will shrink at the rate which is a bit worse than the usual n -1 ^ 2 . 


3 Preliminaries 

3.1 Variation 

In this subsection we remind several facts about p-variation. For details we refer the reader 
to [BJ. Recall that a function / : [a, b] —> R has p-variation defined by 


MI; [MD = «ipEI/(*‘) ~ K^- 1 )T’ 

* k =l 

where p € ( 0 , oo) is a fixed number and supremum is taken over all possible partitions 
x : a = xo < X\ < ■ ■ ■ < x n = b, n ^ 1. 


If v p {f\ [a, 6 ]) < oo, / is said to have a bounded p-variation on [o, 6 ]. Denote by W p ([a, b]) the 
class of functions defined on the [a, 6 ] with bounded p-variation. Let CW p ([a,b]) denote the 
class of functions that have bounded p-variation and are continuous. Define V p (f',[a,b]) = 
v p p (f; [a, 6 ]), which is a seminorm on W p ([a,b]) provided p ^ 1 and V p (f\ [a, b]) is 0 if and 
only if / is a constant. 

If / € W p ([a, b]), then f is bounded and / € W P1 ([a, 6]) for every pi > p > 1. Let 0 < p < oo 
and f,g £ W p ({a,b]). Then fg £ W p ([a, &]). 

Let / £ W q ([a,b\) and h £ W p ([a,bj) with 0 < p < oo, q > 0, 1/p + 1/q > 1. Then 
the integral f b fdh exists as the Riemann-Stieltjes integral if / and h have no common 
discontinuities. If the integral exists, the Love-Young inequality 



dh- f{y)[h(b) - h(a)] 


<C p , q V q (f-,[a,b])V p (h-,[a,b}) 


( 4 ) 


holds for all y £ [a, 6 ], where C p , 9 = £(p 1 + q 4 ) and £(s) = ^J n >i n 
Let / £ Wq([o, b]) and h £ W p ([a, b]). From © it follows that 


V p \ fdh;[a,b] 


< C p , q V qi00 (/; [a, b])V p (h\ [a,b]), 
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where V qt0o (/; [a, b]) = V q (f; [a, b}) + sup^^ |/(x)|. Note that V q , ao (f', [a, b]) is a norm on 
Wq([a,b]), q> 1. 

If h £ CW p ([a, b]), then the indefinite integral f y f dh, y £ [a, b], is a continuous function. 

To verify that some process solves the equation |l| we use a chain rule for the composition 
of a smooth function and a function of bounded p- variation with 0 < p < 2. The chain rule 
is based on the Riemann-Stieltjes integrals. As usual, the composition jo ft on [a, b] of two 
functions g and h is defined by ( g o h)(x) = g(h(x)) whenever h is defined on [a, b] and g on 
the range of h. 

Theorem 4 (Chain rule (see [6])) Let f = (fi, ■ ■ ■, fd) : [a, b] — > R d be a function such 
that, for each k = 1,... ,d, fk S CW p ([a ,b]) with 1 ^ p < 2. Let g : R d —► R be a differentiable 
function with locally Lipschitz partial derivatives g' k , k = 1,..., d. Then each g[of is Riemann- 
Stieltjes integrable with respect to fk and 

( 9 ° f)(b) ~(g° f)(a) = E [ (g'k of) dfk- 


Proposition 5 (Substitution rule (see |B])) Let f, g, and h be functions in CW p ([a,b]), 
1 ^ p < 2. Then 


/ b / nx \ rb 

f(x) d yj g(y) dh (y)j=J 


f(x)g(x)dh(x). 


3.2 Several results on fBm 

Recall that fBm B H = {B : t ^ 0} with Hurst index H G (0,1) is a real-valued continuous 
centered Gaussian process with covariance given by 

T7U tdH r,H\ 1 / 2H . ,2H I, |2 H\ 

E (B t B s ) = -(s +t -\t-s\ ). 

For consideration of strong consistency and asymptotic normality of the given estimators we 
need several facts about B H . 

Limit results. Let 


„ 2 £ T -1 "-1 


Vn.r = 


4 — 2 2H 


E( T ' 






1 

2 ' 


Then (see [3], £)j. [2]) 


with 


o 2 (H) 


V n ,T 

Vi n.T 


- 1 
- 1 


V n ,T - > 1 

n—¥ oo 



a.s.; 

( * 2 (H) 
V ffi.2 (H) 


01 . 2 (H) \\ 
o 2 (H)/2 JJ 


2 + c 2 (H) + ci (H) J2 {p 2 - 2 H (l)) 2 , err, 2 (H) = 2~ 2H (3a 2 + af + dal), 

1=2 


where 


ci(H) 


p2-2H(l ) 


2 


o-i 


( 2H (2H — T)(2H — 2)(2H — 3) \ 2 / 2 2h+2_ 7 _ 3 2h 

4 — 2 2H ) ’ C2 ^ ^ 4- 2 2H 

\l - 2\ 2H - d\l - 1\ 2H + 6\l\ 2H - d\l + 1\ 2H + \l + 2\ 2H ^ 

2H(1 - 2H)(2 - 2H)(3 - 2H) ’ G ’ 

- ° 2 ^ + Cl(H) ^2 P2-2 h(1) p2—2H (l - 2), 


2 


o\ —2\J Cl(H) + Cl(H) ^2 P2—2H (l) P2-2h(1 ~ 1)- 
1=2 
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These relationships imply that 


H a.s., 

Variation of B H . It is known that almost all sample paths of an fBm B H are locally Holder 
of order strictly less than H, 0 < H < 1. To be more precise, for all 0 < e < H and T > 0, 
there exists a nonnegative random variable G £ ,t such that E(|G e ,T| p ) < oo for all p ^ 1, and 

\B? -B?\<G e , T \t-s\ H ~ e a.s. (5) 

for all s,te [0, T], Thus B H e CW He {[ 0,T]), H e = 

SDE driven by fBm B H . A solution of the stochastic integral equation (JT]) , on a given filtered 
probability space (fi, T, P, F) and with respect to the fixed fBm (B H , F), 1/2 < H < 1, and 
initial condition £, is an adapted to the filtration F continuous process X = {Xt : 0 < t < T} 
such that A'o = £ a.s., P(fg \f(X a )\ ds + \ f* g(X a ) dB^\ < oo) = 1 for every 0 < t < T, and 
its almost all sample paths satisfy 0 - 

4 Solution of SDE with polynomial drift 

First, we consider a non-random integral equation 

xt — x o + f (axT + bx s ) ds + c f x s dh 3 , xo > 0, a ^ 0, t € [0, T], ( 6 ) 

Jo Jo 

where to € N, m > 2, h £ CW p ([0, T]), 1 < p < 2. 

4.1 Uniqueness of the solution of equation (|6]) 

Theorem 6 The integral equation 0 has a unique solution in the class of functions GW p ([0, T]), 

1 < p < 2 . 

Proof. It is well-known that 

m — 1 

m i _ra / n \ ^ kim—1 — k 

a — b = (a — b) y. a o 
k =0 

Let us assume that we have two solutions Xt and y t of equation 0 . Since the solutions belong 
to the class of functions GW p ([0, T]), 1 < p < 2, then 

\x? - yT\ ^ mL Xi y trn ,T\xt - yt |, 

where 

U, P ,m,T = (max | max \x t \, max |y t |l) 

V LoOsJt^T J / 

Further, one can find a subdivision r = {0 < ri < T 2 < ■ • ■ < t„ = T} such that 

V p (h-, [Tk-i,T k ]) ^ (4|c|G p , p ) -1 
for all k. Assume we have proved that x Tk _ 1 = Vr k _ 1 - Then 

Vp,oo{x y\ \j~k — 1 ) Dc]) =Vp,oo (x y {x Tk _ 1 2 /rfc_i )? \xk — 1 ? Dc]) 

< 2 |a| f \x?-y?\dt + 2\b\ [ \x t -y t \dt 

dr k _ i 

T 2|c|C'p, p lp,oo(^ y? [De—i? Dc])Tp(/i 5 [r fe _i, Tfc]) 

r T k 

^{2L X: y } m : T\a\ + 2 | 6 |) / |a;t — y t | dt 
J T k -1 

+ 2|c|G p , p V r p ,oo(a; - y; [t*_i, Tfc])tp(/i; [Tfc_i,Tfc]) 


H n = *-— In 

2 2 In 2 


/ Vi n,T 

l22H-iy„, T 


21n2Vn(H n -H r ) ^V(0;a^(H)), cr^tf) = 


5 




and 


Vp ? QO 


V, [Tk-l,Tk]) ^4 (L, 

^4 (L: 


x,y,m,T\Ci\ + 
x,y,m,T |a| + 


r 

1*1) / 

J Tfc 

1 * 1 ) / 

J T 


T k -1 
Tfc 


|x f - 2/t| dt 


Vp ? oo (* - y, [r k -i,t])dt. 


By Gronwall’s inequality we then have V Viao {x — y\ [ 7 >_i,Tfc]) = 0. Thus x t = yt on [r k - 1 ,Tk\- 
Since x TQ = xo — y T0 the claim of the theorem follows from repetitive application of the above 
reasoning. 


4.2 Explicit solution of equation (|6]) 

Theorem 7 Function 

/ ft \ i/l 1 -™) 

Xt = e u+cht (xl- m + {l-m)aJ e ( — 1)( ^+^) ds \ , te[0,T], (7) 

is an element of the class of functions CW p ([0, T]), 1 < p < 2, and satisfies equation 0 - 

Proof. We first show that x G CW p ([0, T]), 1 < p < 2. Let us denote Zt = exp {bt + cht } and 
z k = Zt k , where {0 = to < ti < ■ • ■ < t n = T} is a subdivision of the interval [0,T]. By the 
mean value theorem for the function e x , x G R, we obtain 

n —1 n —1 

X) \ Zk + x ~ Zk \ P ^ ex P{p[l & l' r + 2 |c|V p (/i; [ 0 ,T])]} l bA fc+i + cAhk+i\ P 

k =0 k =0 

^2 exp {p[| 6 |T + 2|c| Vj, (h; [0, T])] } [| 6 | P T P + |c| p u p (/i; [0, T])], 
where A fe+i = t k+1 - t k and Ah k+1 = ht k+1 - ht k . Let 

ft = xl~ m + (1 - m)a [' e (m - 1 )(i,s+cM ds. 

Jo 

Repeatedly using the mean value theorem for the function u H>■ u u > 0, m ^ 2, we 
obtain 


n — 1 

k =0 


/ fc 1/(1 " m) l £ IA+i - AI < *0^1 (/; [o,n)- 

fc =0 


Since z G CW p ([0,T]) and / G CWi([0,T]), we have x G C'W’ p ([0,7 1 ]). 

Now we verify that function 0 satisfies the equation (j6j by usingthe chain rule. Indeed, 
let F(t,u,v) = e bt+cu v 1/ '^ 1 ~ m K Then by Theorem [ 4 ] and Proposition [5] we get 


x t =F(t, h t , ft) = F( 0, h 0 , fo) + b / x s ds + c / x s dh. s 
1 rt 


if Xsds + c 

Jo Jo 


1 — m 


/ 


e bs+ch s j l/(l-m)-l d ^ 


=xo + b Xsds + c x s dh s + —-— 

Jo Jo 1 - m J 0 


bs-\-ch s rTn /(1 —m) 

Js 


dfs 


=xo + b f x s ds + c f x s dh s + -—-— f (1 — m)ae rn ^ bs+ch ^ f™^ 1 
Jo Jo 1 - rn J 0 


ds 


i Xsds + c 

Jo Jo 


=xo + b Xsds + c x s dh s + a / x™ ds. 


4.3 Solution of SDE 

Since almost all sample paths of B H , 1/2 < H < 1, are continuous and have bounded Ft e = 
-g£-variation, 1/2 < Ft — e < 1, the pathwise Riemann-Stieltjes integral J* X s dB f exists for 
X G CWjf E ([0,T]). So SDE 0 makes sense for almost all uj. Thus the obtained result for a 
non-random integral equation can be applied to an equation driven by a fBm. 
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Theorem 8 Suppose that Xq > 0 and m ^ 2. The stochastic process 

V- bt+cB¥ f \rl — m , /-I \ f (m — l)(bs-\-cBT) j 

X t = e * I X 0 + (1 — m)a e K s J ds 

for almost all uj belongs to CWh £ ([ 0,T]) and is the unique solution of 0 - 


\ l/(!-m) 


t£[0,T\, 


(8) 


5 Proof of the main Theorems 


5.1 Proof of Theorem [Q 

Before presenting the proof of this theorem, we give an auxiliary lemma. To avoid cumbersome 
expressions we extend notions of o r , O r to continuous time processes as follows: Yt = O r (q(t)) 
means that there exists some almost surely finite non-negative r.v. £ and q(t) -> 0 +, with 

the property \Yt\ (q(t ), t — > 0, whereas Yt = o r {q(t)) means that \Yt\ ^ £«;(£), t —¥ 0, with 
w(t ) = o(q(t)), t —> 0 . 

For t £ [0, T ) and li^Owe also denote 

A Y t ,t+ h = Y t+h - Y t , A (2) Y t ,t+h = A Y t+h}t+2h - A Y t , t+h . 

Lemma 9 Assume that H £ (1/2,1), a ^ 0, m ^ 2 and X is the solution of equation 0 - 
TTien for every e > 0 such that H > 1/2 + e 

&Xt,t+h = Xt(bh + cAB t ,t+h + O r (Ji )), h — > 0+, t£[0,T), 

and 

A^Xt't+h = X t (cA {2) B B t+h + O r (h 2(H - e) )), h -»■ 0+, t £ [0,T). 

Proof. We rewrite the solution of equation (jT]) in the following form 

X* = A' 0 Z t A t , 


where 

Zt = exp{ 6 £ + cB b } and = ^1 + (1 — m)aX™~ 1 J Z™~ 1 ds ^ 
Let h —£ 0+. By ([5| and Maclaurin’s expansion of e x 

Z t+h me Kt+h) e cB ?+ h = e b \l + bh + o(h 2(H ~ e)) )e cB ^ +cAB ^ 
=Z t (l + bh + o(h 2{H ~ e) ) s j (l + cAB^+h + O r (h 2(H ~ s> )) 
=Z t (l + bh + cAB B t+h + O r (h 2< - H ~ e) )) . 

Next note that t t—> At ^ 1 is non-decreasing whereas 

rt-\-h 


l/(l-m) 


/ 


l ds = hZ™- e f l =Or{h), 6 £ (0; 1), 


vanishes. Therefore Maclaurin’s expansion of x H ¥ (1 + x) a gives 

r t+h \ 1 /( 1 - m ) 


/ B 

A t+ h = P4 t 1-m + (1 - m)aX™ _1 J 

=At (l + It,t+h + Or(h 2(H ~ e) )) , 


Zf l ~ 1 ds 


where 


Hence 




aX, 

~A 


rm — 1 rt-\-h 

1 — m I 
t Jt 


X t+h =X t (l + bh + cAB B t+h + O r (h 2(w - £) )) (l + I tit+h + o r (h 2(H ~ s) )) 
—Xt (l + bh + cAB B +h + It,t+ h + 0.(h 2( "“ £) )) . 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


7 




Consequently, 


AX t ,t+h = X t (bh + cAB” +h + I t , t+h + O r (h 2(H ~ e) )) 

and taking into account (121 it follows that 

A {2) X t ,t+h =X t+h (bh + cAB? +h ' t+2h + I t+h ,t+2h + O r {h 2(H ~ s) )) 

- X t (bh + cAB* +h + I t , t+h + Or(h 2(H ~ s) )) 

=A't(l + bh + cAB^ t+h + It,t+h + O r {h 2 ^ H e ')) 

X (bh + cAB t +h,t+ 2 h + It+h,t+ 2 h + Or{h * E '*)) 

- Xt (bh + cAB* t+h + It,t+h + O r (h 2{H ~ E) )) 

=Xt (cA (2) B» t+h + I t +h,t+2h - It,t+h + Or(h 2{H - £ >)) . 

It remains to show that I t +h t+ 2 h — It t+h = O r (h 2 ^ H ~^). To see this observe that by 
@ (JTTj) and the mean value theorem, for some 6t £ (0; 1), * = 1,2, 


It+h,t+2h It,t+h — 


aX, 

A 


m — 1 / rt-\-2h 

izzr ( / 

+ h \Jt- 

UU 

;+h \J t- 


-1 rt+h 


i t+h \Jt + h 

aX r m_1 f r t+2h 


1 t+h \Jt+h 

a xr ~i / rt+2h 


Tl( [ 

r+h \Jt+h 


z\ 

z m ~ 


•-mi 

1 ds-(l+Or{h )) rn ~ 1 j 


Z™~ 1 ds + O r (h 2 ) 


aXZ 


A 


- mm+eth - hZr+elh + O r (h 2 )) 


t-\-h 


aXZ 


A 


t-\-h 


hz: 


m — l I I ^s+t 


- 1 +O r (h 2 ) 


aX, 


(hZT - 1 ((l +O r (h H - e )) m 1 - l) +Or(h 2 )) 


=^ro r mm=o r m H -% 

A t+h 

where s = t + d 2 h and t = h( 1 + 9i — 0 2 ). 

Proof of Theorem [7j Observe first that the function 

/ 'J 1 \2x 


<{>n,T(x) = ( —) (4 — 2 2X ), X £ (0, 1), 


is continuous and strictly decreasing for n > T. Thus it has the inverse function ip n 'r f° r 
n > T. By Lemma [9] we obtain 


tpn,T(H1 






Tti,t(H) 


=} / a (2) 


cX k 




- 1 / A (2) 


cX k 


T 2H (4- 2 2H ) 


~2H-1 


n 

~4-2 2H 


EC 


X k= 1 7 

y+o,^),v n , T+ o,.(^ rs ) 


t~ h a^ ) b ' 


for 3e < H. 
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Fix 5 > 0. Let /3 n = ■ Since V n ,t 1 a.s., the same applies to /3 n . Consequently, 


Vn,rW 

there exists ng(uj ) such that 


<Pu,t(H + 5) =(4 - 2 2(j?+5) ) Q 


y\ 2 (fl+S) 


</3n(4-2' J ") ( - 

' n 


/ r r\ 2(H — 5) 
<(4-2 2(h - s) )^-) = Vn, T (H — 8) 


for all n ^ no, i.e. Pu,t(H^) £ (<pn,T(H + 8),<pn,T{H — 5 )) for all n ^ no- Since the 
function <p n ,T is strictly decreasing for all n > T we conclude that H„ 1> £ ( H — S, H + <5) for 
all n ^ max(no,T). Thus H © is strongly consistent. 

Now we prove the asymptotic normality of H^\ Note that 


In 


(i)> 


P>n,r(H) 


= In ( — 


© 


-2 


+ In 


4-2 


Hi) 


4-2 2H 


= - 2(J?i 1) - H) In (p) + In (4 - 2 2H - In (4 - 2 2 "). (13) 

By the Lagrange theorem for the function h(x) = ln(4 — 2 2x ), x £ (0,1), we obtain 
In (4 - 2 2S " 1> ) - In (4 - 2 2H ) = ( Hf ) - H)ti (H + 9 n (H « - H)) 
for some 9„ £ (0,1). Since Hn' 1 —> H the sequence ( h’(H + 9 n (H ,© — H))) n j,i is bounded. 


The equality (131 can be rewritten in the following way: 

( 


-2hi (£«-#) 


n( In I 


v n , T (H^) 

<Pn,T( H ) 




h'(fr+e rl (S( l 1, - 
21n (t) 


n ( ln ( ^Cm ) )- lnl ) 

1 + O r (l) 


(14) 


Next note that the Delta method together with Slutsky’s theorem and limit results of section 
|3.2| imply 




l fin,T(H i 

P>ti,t(H) 


l) = Vs(v„,t - l) + Or( nfl _3._ 1/3 ) N(0-,a 2 (H)) 


once £ > 0 is chosen to satisfy 3e < H — 1/2. It remains to observe that assertion |2]| imme¬ 
diately follows from (141 and repeated application of the Delta method along with Slutsky’s 
theorem. 


5.2 Proof of Theorem [ 

Fix e £ ^0, • By Lemmajoj 


E 


Xn.k 


-J2(cA %b» + o 


,-2(H-e)\y 


k=1 
n — 1 


; E( A i>T+° 


1-3 (H-e) 


k=1 
J2rr2H 


2r-p2H (a r\2M \ 

i© T .., r+0 . („-»■«-») 


and in the same way 

2n-l / a (2 ) 


E 


Xnn.k 


c 2 T 2H ( 4-2 2H ) 
(2n) 2H ~ 1 


V2 n ,T + O r (n 1 - 3{H - e) ). 
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Since V„,T,V 2 n,T —> 1 a.s. (see section |3.2 1 
1 


' 2 |M W+ 0 r(n .-3(H- £ ))' 


H, 


( 2 ) 


In 


2 2 In 2 I ^( 4 _- r) KT + a(nl — 3(g — £)) 


1 


In 


H 2 n,T / 1 + O r (n _W+3e ) 


2 2 In 2 1 2 2if - 1 K, T l l + O r (n-»+3E) 


—~B n 


2 In 2 


In (l + O r (n- w+3e )) = H n + O r (n~ H+Z£ ). 


Now to finish the proof apply Slutsky’s theorern_and limit results of section [372| Note that the 
limit variance cr 2 (if) of H „ equals to that of H n . 

5.3 Proof of Theorem |3] 

Let 0 < £ < 1/2 be such that 3 e < H — 1/2 and 

^ 2 n 2 ("- S " 3 ) )( 4 - 2 2S " 3) ))^ 2 

On * fQi 0 n . 

_ 2 2-ff) 

By Maclaurin’s expansion and the mean value theorem we obtain 


1 2 ( m ( 3 >)( 4 2 2 »' 3) )) 


J-2(H-h/ 3) )(4 _ 2 2H) 


= exp |o r ((/(n)) In + In (4 - 2 2W " )-ln(4-2 2H )| 

=(1 + o r {<t>{n) Inn)) exp {o r {<j>(n))h! (H + 0„(H^ — H))} 
=(1 + o r (/>(n) Inn))(l + o r ((/>{n))) = (1 + o r (</(n) Inn)) 


for some 8 n £ (0,1), where h(x) = ln(4 — 2 2x ), x £ (0,1). Hence it suffices to show that c/ 
satisfies 

c n —> c a.s., v n(c n — c ) —> 7V(0; c a ) 
and the claim will be proved. By Lemma [9] 

2 n — 1 


E 




Afe A ' -'’ 2 E( A i 2)B ") 2 + a ( nl " 3(H " e) )- 


Thus 


2 2 / 
— r —c. 


2H-1 n-1 


n — i \ 

Y (A (2) sf) 2 - 1) +O r (n- H+3e ) 
fc = l ' 


VT 2K (4-2 2i? ) 

=c\V n , T - 1) +O r (n~ H+3s ). 

It remains to apply Slutsky’s theorem and limit results stated in section [372] 


6 Simulations 

The simulations of the obtained estimators presented below were performed using the R soft¬ 
ware environment m- The performance of these estimators was assessed using sample paths 
of the Black-Scholes model, the Verhulst equation and the Landau-Ginzburg equation as de¬ 
fined below. In Equation [I] if we set a = 0, b = A, c = a we obtain the Black-Scholes 
model 

dX t = XX t dt + aX t dB f. 

Let us note that the results, obtained in Section [2] are applicable for the Black-Scholes model. 
If we set a = —1, b = \, c = a,m = 2 we obtain the Verhulst equation 

dX t =(- X 2 + A X t )dt + aX t dB f 

with the explicit solution 

x = gexp{A t + o-Rf } 

1 + £ exp{As + aBfjds 
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If we set a = —1, b = X, c = a, m = 3 we obtain the Landau-Ginzburg equation 


dX t = ( - X? + XX t )dt + aX t dB f. 


with the explicit solution 


£ exp{At + crBf} 

1 + 2£ 2 f* exp{2(As + aBp)}ds 

Further on, unless explicitly stated otherwise, the (arbitrary) default values of £ = 3, 
b = 0.5 and c = 0.7 were used, and in each case the conclusions were drawn from 500 sample 
paths of the respective processes. The figures discussed further can be found in the Appendix. 

Firstly in Figure[I]we consider the dependance of numeric characteristics of the estimators 
H { ,y (defined in Theorem [l] and denoted in the figures as HI) and Hi 2> (defined in Theorem 

[2] and denoted in the figures as H2) on the true value of the parameter H. 

The results imply that the accuracy of surpasses that of by at least an order of 
magnitude, and becomes even more precise for the values of H close to 1. The practical use of 
Hn\ however, requires the knowledge of the true value of the parameter c, which is not always 
available. However if that is the case then the estimator Hit' yields accurate estimators of H 
even for very short sample paths, as illustrated in Figure [2] It is also worth to note that the 
estimators do not display notable dependance on the underlying process. 

The next question on the agenda is, how rapidly does the accuracy of the estimators grow 
as the length of sample paths is increased? Figure [3] presents the boxplots of both estimators 
as the length of the sample paths, N, varies from 1024 to 8192 points. It can be seen that the 
estimator is considerably more precise than H„\ yet the rate at which their accuracy 
grows with increased sample sizes is roughly the same. For both of these estimators the 
inter-quartile ranges shrink by roughly 30% as the sample length is doubled. 

Now let’s examine the behavior of ? 2 , the estimator of the volatility defined in Theorem [ 3 ] 
and denoted c2 in the figures. As before, we present the boxplots of this estimator both for 
different values of the true parameter c (Figure |d| and different sample path lengths (Figure 
[5j. We can observe that the distribution of the estimator is somewhat right-skewed, which is 
hardly surprising since we’re estimating c 2 , not c itself. Also, Figure [ 4 ] shows that the variance 
of c 2 grows rapidly as the true value of the parameter increases. This is further illustrated by 
Table nl presenting the average bias S' 2 — c 2 along with the average variance of this estimator 
for different values of c. The adjusted R-squared of the linear model Var(ct ) = kc 4 + b is 
0.9838 which is in agreement with the asymptotics presented in Theorem [3] 



c 

0.2 

0.5 

1 

2 

5 

Bias 

0.006 

0.031 

0.133 

0.460 

3.188 

Variance 

0.001 

0.019 

0.306 

5.159 

193.4 


Table 1: Numeric characteristics of 
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Appendix: Figures 
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Figure 1: The differences between the estimated and real values of H 
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Figure 3: The differences between the estimated and real values of H 
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Figure 4: c„ — c 2 , dependance on c 
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